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ABSTRACT 

Recent observations on the cosmic microwave background by Wilkinson Microwave Anisotropy Probe 
(WMAP) strongly suggest that the reionization of the universe took place quite early {z ~ 17). On the 
other hand, it has been pointed out that the cold dark matter cosmology suffers from the substructure 
problem that subgalactic halos are overproduced than the observed dwarf galaxies in the Local Group. 
In this paper, as a potential mechanism to solve this problem, we consider the feedback effects of the 
early reionization on the formation of small-scale structures. For this purpose, we perform 3D radiation 
hydrodynamic simulations with incorporating the radiative transfer for ionizing photons. As a result, 
it is found that the early reionization is so devastating for low-mass systems with M^h- < IO^Mq or 
''^circ ^ 20km s~^, and almost all gas is photo-evaporated in more than 95% of low-mass systems. Such 
a strong negative feedback on the formation of low-mass galaxies may solve the substructure problem 
and support the picture that Local Group dwarf galaxies are descendants of the more massive halos that 
experienced and survived tidal stripping . 

Subject headings: galaxies: formation — galaxies: dwarf — radiative transfer — molecular processes — 
hydrodynamics 



1. INTRODUCTION 

It has been claimed that too many subgalactic halos are 
formed in a cold dark matter (CDM) universe, compared 
to the dwarf galaxies observed in the Local Group (Klypin 
et al. 1999; Moore et al. 1999). One of the simplest ways to 
resolve this discrepancy is the exclusion of baryonic mat- 
ter from small halos by some feedback mechanisms. The 
mechanisms could be either the gas ejection due to energy 
input by the multiple supernova (SN) explosions (Dekel & 
Silk 1986; Yepes et al. 1997; Efstathiou 2000; Kay et al. 
2002; Marri & White 2002; Wada & Venkatesan 2003; Ri- 
cotti & Ostriker 2004) or the photo-evaporation by ultra- 
violet background (UVB) radiation (Umemura & Ikeuchi 
1984; Bond, Szalay & Silk 1988; Efstathiou 1992; Babul 
& Rees 1992; Thoul & Weinberg 1996; Barkana & Loeb 
1999; Kitayama & Ikeuchi 2000; Kitayama et al. 2000, 
2001; Bullock et al. 2000; Somerville 2002; Benson et al. 
2002; Ricotti, Gnedin, & Shull 2002; Susa & Umemura 
2004). As for the former process, Wada & Venkatesan 
(2003) recently performed high resolution hydrodynamic 
simulations and concluded that 1000 SNe are required in 
order to disrupt a galaxy with M ~ IO^Mq at z ~ 10, 
while only 100 SNe in the halo lead to triggering the re- 
collapse of the whole system. These results show that SN 
explosions are not always destructive for the formation of 
low-mass galaxies, although it depends upon the initial 
mass function of stars (e.g. Ricotti & Ostriker 2004). On 
the other hand, the previous studies on the UVB feedback 
have shown that the photo-evaporation by UVB is depen- 



dent upon the virial mass. As long as systems are more 
massive than w 1O^M0, the photo-evaporation by UVB is 
not effective, because deep potential wells can retain the 
ionized gas. However, at smaller mass-scales, the feedback 
by UVB is expected to play a very important role. Re- 
cently, Wilkinson Microwave Anisotropy Probe ( WMAP) 
(Kogut et al. 2003) suggested that the universe was reion- 
ized in a rather early epoch (z ~ 17). If this is the case, 
the UVB feedback could be quite effective for low-mass 
systems, since density fluctuations are photo-heated be- 
fore they collapse to form stars. 

Very recently, Kravtsov, Gnedin & Klypin (2004) an- 
alyzed the dynamical history of small substructure halos 
with present mass < 10^~^Mq by a cosmological simula- 
tion. They found that 10% of small halos originate in con- 
siderably larger systems with > IO^Mq, which survived 
the tidal stripping. They suggest that the Galactic satel- 
lites are descendants of relatively massive systems which 
formed at higher redshifts and did not significantly suf- 
fer from the UVB. This argument is plausible, if systems 
with > IO^Mq are impervious even to the early reioniza- 
tion, and also the star formation in low-mass systems is 
completely suppressed. 

Kitayama et al. (2001) have investigated the UVB feed- 
back by one dimensional radiation hydrodynamic simula- 
tion for halos which collapse at z < 10, and found that the 
central parts cannot cool down to form stars if the virial 
temperature is less than 10* K and I21 > 10~^, where 
/21 is the intensity at Lyman limit in units of 10~^^ erg 
str~^. Dijkstra et al. (2004) extended this 



work to higher redshifts by one dimensional hydrodynamic 
calculations, and found that gas is not photo-evaporated 
in halos with circular velocities of 10 — 20km s~^ at z > 10. 
But, in the hierarchical structure formation in a cold dark 
matter universe, three-dimensional (3D) time-dependent 
self-shielding is significant for the UVB feedback, as shown 
by Susa & Umemura (2004). In particular, such radiative 
transfer effects are likely to be essential, when the rcion- 
ization proceeds at higher redshifts (in higher density in- 
tergalactic matter). 

In this paper, we perform 3D radiation hydrodynamic 
simulations with solving the radiative transfer for ioniz- 
ing photons, and attempt to elucidate the feedback effects 
of the early reionization on the formation of subgalactic 
systems. In the next section, the method of numerical 
simulations is briefly described, and the assumption for 
ionizing radiation intensity is provided. In §3, we present 
the numerical results for UVB feedback and also make the 
convergence check of the results. §4 is devoted to the dis- 
cussion on the substructure problem. 

2. NUMERICAL METHOD 

The details of numerical method is given in Susa & 
Umemura (2004). Here, we briefly describe the model and 
method for simulations. 

Hydrodynamics is calculated by Smoothed Particle Hy- 
drodynamics (SPH) method. We use the version of SPH by 
Umemura (1993) with the modification according to Stein- 
metz & Miiller (1993), and also we adopt the particle resiz- 
ing formalism by Thacker et al. (2000). The gravitational 
force is calculated by a special purpose processor for grav- 
ity, GRAPE-6 (Makino 2002). In order to access GRAPE 
boards, we utilize the Heterogeneous Multi-Computer Sys- 
tem (HMCS) (Boku et al. 2002) which allows us to use 
GRAPE in parallel processors such as PC clusters. 

The softening length for gravity is set to be 20pc for SPH 
and CDM particles. This value is based on the convergence 
test shown in section 3. The non-equilibrium chemistry 
and radiative cooling for primordial gas are included with 
the code developed by Susa & Kitayama (2000), where 
H2 cooling and reaction rates are taken from Galli & Palla 
(1998). To solve radiative transfer in an SPH scheme, we 
employ the method proposed by Kessel-Deynet & Burkert 
(2000), which utilizes the neighbour lists of SPH particles 
to evaluate the optical depth from a certain source to an 
SPH particle. 

To generate UV background radiation, we put a single 
ionizing source located very far from the simulated region, 
and control the UV intensity by specifying the incident 
flux to the simulation box. We assume the history of 
UV intensity that allows the early reionization inferred by 
WMAP. Since the recombination time scale is shorter than 
the Hubble expansion time at 2: > 10 (Peebles 1968), ion- 
izing photons should be continuously supplied in order to 
retain the ionization. The requisite UV intensity is given 
by equating the recombination rate to ionization rate: 

no(l + zfk,echp{a + 3) (1 - ym)^ 



Iu^ = 



47rcr„ 



VBl 



(1) 



where Ii,^ denotes the intensity at Lyman limit, and fcrec, 
cr,yL 1 ^iid ^P respectively represent the recombination coef- 
ficient, ionization cross-section at Lyman limit, and Plank 
constant, no is the present-day baryon density, a is the 



spectral index of UV radiation, and j/hi is the fraction of 
neutral hydrogen. If yni = 0.1 at z ~ 15, I21 — 10^*^ is 
required. This should be regarded as the minimal prereq- 
uisite for two reasons. First, this estimate is provided for 
optically-thin media. In actual ionization history, the ra- 
diative transfer effect is definitely important (Nakamoto, 
Umemura, & Susa 2001). Secondly, in an inhomogeneous 
universe, local density enhancements increase the recom- 
bination rate significantly. Taking these into account, we 
assume I21 ~ 0.01 for 5 < z < 17 with a sharp decrease at 
z > 17 as /21 oc exp [3 (17 — z)]. Also, some models with 
/21 = 1 are simulated to investigate the dependence on 

hi- 

The "star formation" condition adopted in this paper 
is basically the same as Susa & Umemura (2004), except 
that c* = 1 is assumed here, c* = 1 means the star forma- 
tion at maximal feasible rate, because no stars can form 
in a timescale shorter than local free-fall time. Here, we 
do not include the internal feedback effects by stellar UV 
radiation as well as SN explosions. 

We assume ^m = 0.3, f^A = 0.7, f^baryon/i^ = 0.02, and 
h = 0.7, as the background cosmology. Following the den- 
sity fluctuations in this cosmology, the initial distributions 
of particles are set up. The mass of virialized dark halo is 
in the range of IQ^Mq < M^-„ < IO^Mq, and the coUapse 
epoch is 5 < Zc < 20. We use 2^^ SPH particles and 2^^ 
dark matter particles for a run. 

3. RESULTS 

The numerical results are summarized in Figure 1. The 
left panel shows the fraction of the final mass in stellar 
component to the initial baryon mass as a function of the 
collapse epoch z^ and Mvir, while the right panel is the 
final fraction of gas. In both panels, the fraction is de- 
picted by different symbols. Dotted lines represent the 
collapse epoch of halos formed from Icr, 2a and 3ct CDM 
density fluctuations. Roughly 95% of fluctuations collapse 
after the epoch predicted by the 2cr line. Dashed lines de- 
note the constant circular velocities, which arc defined by 
Wcirc = \/GMvii/rvir with the virial radius Tvir determined 
by Zc and Mvir- 

In Figure 1, we see that, if Mvir ^ IO^Mq and Wcirc }t 
20km s~^, almost all baryons are transformed into stars, 
leaving little gas. It is noted that a considerable fraction 
of baryonic matter form stars even after the reionization 
(zc < Zreion = 17). On the other hand, if Wcirc ^ 20km s^^, 
/star steeply decreases with decreasing circular velocities. 
Also, /gas becomes quite small as seen in the right panel. 
This reflects the effect of photo-evaporation by the UVB. 
The dependence of photo-evaporation on Mvir and z^ is 
understood by the self-shielding and the gravitational po- 
tential. The self-shielding against UVB becomes promi- 
nent if the local density exceeds a threshold density as 
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(Tajiri & Umemura 1998). This means that the self- 
shielding is more effective if the system is larger or col- 
lapses at higher redshifts. The gas envelope that is not 
shielded from UVB is photo-heated to < 10'*K and is blown 
out by the enhanced thermal pressure, unless the gravita- 
tional potential is deep enough to retain the ionized gas 



(i.e. iicirc ^ 20kni s~^). Hence, the photo-evaporation is 
quite devastating for fluctuations with lower masses and 
later collapse epochs. Intriguingly, a steep transition of 
/star for Wcirc ^ 20kms~^ coincidentally lies on a line of 
nearly constant cr, i.e., 5p/ p « 2a. This means that in 
more than 95% of the halos with Wcirc ^ 20kms^^, the star 
formation is strongly suppressed by the early reionization. 

In the present simulation, the star formation is assumed 
to proceed in local free-fall time. This leads to the phys- 
ically maximal star formation rate. Hence, the obtained 
stellar fraction in a halo should be regarded as a maximal 
one. If the star formation proceeds in a longer timescale, 
baryon gas in a halo with Wcirc ^ 20kms^^ could be al- 
most completely photo-evaporated after the reionization. 
On the other hand, fluctuations with Wcirc <^ 20kms^^ is 
like to be impervious to the photo-evaporation, since the 
gravitational potential is deep enough to retain the ionized 
gas. 

To examine the the dependence of results on assumed 
UVB intensity {Iii), we also perform several runs with 
/21 = 1. In this case, a steep transition of /star is slightly 
shifted to higher a as 5p/ p k, 2.5(t. Thus, it turns out 
that the results are not strongly dependent on the UVB 
intensity. Such insensitiveness may be understood by the 
weak dependence of the self-shielding on the UV intensity 
seen in equation (2). 

Finally, to check the numerical effects, we analyze the 
convergence of the runs with changing the mass resolu- 
tion and softening length. In Figure 2, the left panel 
shows the final stellar fraction (/star) in runs with Mvir = 
6 X 10^ Mq and Zc ^ 10, while the right panel is with 
Mvir = 6 X IO^Mq and Zc ~ 15. They are at the transition 
region of /star in Figure 1. The horizontal axis in Figure 2 
is the softening length e for gravity. Four different curves 
correspond to the different numbers of particles used in 
the simulations. We can see that /star almost converges, if 
-^SPH ^ 2^^ and if e < 10 — 20pc. Thus, the present simu- 
lation with e = 20pc and A^sph = 2^'''(= A^dm) is unlikely 
to suffer from numerical effects. 

4. DISCUSSION 

The present numerical simulations predict strong neg- 
ative feedback on the formation of dwarf galaxies with 
^circ ^ 20km s^^. Based on /star, we can make a rough 



estimation of the mass-to-light ratios (M/L) for finally 
formed galaxies. If we assume Mstar/L = 3 for stars in 
solar units, /star = 0.01 corresponds to M^ir/L — 2.6 x 10^ 
and /star — 0.1 to M^ir/L — 2.6 x 10^. But, these values 
cannot be compared directly with the observed M/L of 
satellite galaxies, because a quite large fraction (typically 
more than 90 percent) of dark matter halos of satellites 
can be tidally stripped, as shown by Kravtsov, Gnedin & 
Klypin (2004). Thus, the eventual M/L of satelhtes is 
likely to decrease by a factor of 10, e.g. M/L — 2.6 x 10^ 
for ./star = 0.01. Local Group dwarf galaxies (dSphs and 
dirrs) exhibit a wide range of M/L, which is from a few 
up to w 100 (van den Bergh 1999; Mateo 1998; Hirashita, 
Takeuchi, & Tamura 1998). Hence, the formed galax- 
ies with /star ^ 0.01 may account for the Local Group 
dwarf galaxies. The present simulation predicts that only 
a few percent of fluctuations result in /star ^ 0.01, if 
Wcirc ^ 20kms~^. If the star formation rate is lower, the 
probability is reduced further. Hence, such intrinsically 
low-mass halos may be too few to account for all Galactic 
satellites. 

One possibility to reconcile this discrepancy could be a 
model suggested by Kravtsov, Gnedin & Klypin (2004). 
They found that 10% of small halos originate in consid- 
erably larger systems with > 1O^M0, which survived the 
tidal stripping, and suggest that the Galactic satellites are 
descendants of relatively massive systems which formed at 
higher redshifts. In our simulation, systems larger than 
1O*M0 are not subject to the UVB feedback. Thus, their 
model seems viable to account for the number of Local 
Group dwarf galaxies. 
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Fig. 1. — Summary of numerical runs are shown as a function of the collapse epoch Zc and the virialized halo mass Afvir- The left panel 
shows the fraction of the final mass in stellar component to the initial baryon mass, while the right panel is the final fraction of gas. In both 
panels, the fraction / is depicted by different symbols; filled circles are / > 0.9, filled squares 0.5 < / < 0.9, filled triangles 0.1 < / < 0.5, 
open triangles 0.01 < / < 0.1, open squares 10~^ < / < 0.01, and open circles f < 10"-^. Dotted lines represent the collapse epoch of halos 
formed from Icr, 2a and Sit CDM density fluctuations. Dashed lines denote the constant circular velocities of 5 km/s, lOkm/s, and 20km/s. 
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Fig. 2. — Convergence check of the numerical simulations. The left panel shows the convergence of runs with M^ 



and the right panel is for M^ 
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6 X 10 MQandzc ~ 15. In both panels, vertical axes represent the final stellar mass fraction, and the 



horizontal axes show the softening length of the gravitational force. Four different curves represent the results with various number of 
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crosses A^SPH = 2 , triangles Ngpu 



and squares Nsph = 2 



